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Abstract 

Lymphotoxin-alpha (LTA) is a pro-inflammatory cytokine that plays an important role in the inflammatory and 
immunologic response. Numerous studies have shown LTA polymorphisms as risk factors for cancers, but the results 
remain inconclusive. The goal of the present meta-analyses is to establish the associations between cancers and 
four LTA variants (rs1041981, rs2239704, rs2229094 and rs 746868). A total of 30 case-control studies involving 
58,649 participants were included in the current meta-analyses. Our results showed significant associations with 
increased cancer risk for rs1041981 (odd ratio (OR) =1.15, 99% confidential interval (CI) = 1.07-1.25, P < 0.0001, I 2 
= 12.2%), rs2239704 (OR = 1.08, 99% CI = 1.01-1.16, P = 0.021, I 2 = 0.0%) and rs2229094 (OR = 1.28, 99% CI = 
1.09-1.50, P = 0.003, I 2 = 0.0%). No evidence was found for the association between rs746868 and cancer risk (OR 
= 1.01, 99% CI = 0.93-1.10, P = 0.771, I 2 = 0.0%). Subgroup meta-analysis suggested that rs2239704 was likely to 
increase the risk of hematological malignancy (OR = 1.10, 99% CI = 1.01-1.20, P = 0.023, I 2 = 0.0%), and rs2229094 
was specific for the increased risk of adenocarcinoma (OR = 1.33, 99% CI = 1.11-1.59, P = 0.002, I 2 = 0.0%). In 
conclusion, our meta-analyses suggested that the LTA rs1041981, rs2239704 and rs2229094 polymorphisms 
contributed to the increased risk of cancers. Future functional studies were needed to clarify the mechanistic roles of 
the three variants in the cancer risk. 
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Introduction 

With the high prevalence and mortality rate, cancers have 
become one of the main causes of morbidity and mortality 
worldwide [1]. There are over 200 different known cancers in 
humans, and the mechanisms of cancer pathogenesis remain 
obscure [2]. Large epidemiologic and clinical studies have 
illustrated that inflammation may be associated with the 
development of cancers [3,4]. Increasing evidence suggests 
that predisposition to cancer is associated with cytokines [5,6], 
such as tumor necrosis factors (TNF) [7,8]. 



Lymphotoxin-alpha (LTA) is a pro-inflammatory cytokine 
belong to the TNF family which plays an important role in the 
inflammatory and immunologic response [9]. LTA is a product 
of stimulated T cells [10,11], and it can help communicate 
lymphocytes and stromal cells and subsequently eliciting 
cytotoxic effects on cancer cells [9,12]. LTA gene is located on 
the 6p21.3 that harbors the class III region of the major 
histocompatibility complex (MHC) locus [13]. Genetic variations 
of inflammation-related genes are shown to alter both the 
regulation of the inflammatory response and modulation of 
susceptibility to radiation-induced normal tissue damage [14]. 
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LTA gene polymorphisms are shown to be associated with the 
inflammatory and immunomodulatory diseases including 
cancers [15,16]. The association is significant between LTA 
polymorphisms and cancers including gastric [15] and breast 
[17] cancers in Asians, colorectal cancer [18] in Germans, and 
Non-Hodgkin Lymphoma (NHL) in Europeans [19-21]. 
However, discrepancies remain for the association of LTA 
polymorphisms with cancers in different ethnic groups 
[15,16,22-24]. 

Four LTA polymorphisms (rs1041981:Thr26Asn; rs2239704; 
rs2229094: Cys13Arg; rs746868) have been extensively 
investigated as potential risk factor for cancer. These mutations 
may exert possible regulatory regulation on the gene 
expression and the level of cytokine production [25,26]. Single 
nucleic polymorphism (SNP) rs1041981 (Thr26Asn) [27] is 
associated with the transcriptional regulation of LTA, which 
may activate the lymphocytes and induce apoptosis [27]. SNP 
rs2239704 can modulate both LTA levels and subsequent 
inflammatory response to pathogens [28]. SNP rs746868 is 
shown to be in high linkage disequilibrium (LD) with rs1041981 
(r 2 > 0.8) [29]. SNPs rs746868 and rs2239704 are significant 
predictive variables of LTA protein production [28]. Haplotype 
of SNP rs2229094 is shown to be associated with altered LTA 
expression and increased levels of vascular- and autoimmune- 
mediated inflammation [30]. 

The associations of LTA variations with cancer risk have 
been evaluated by several case-control studies [18-20,31-52]. 
Previous studies have suggested that the functional 
polymorphism rs909253 is associated with gastric [15] and 
breast [17] cancers in Asians. Although the four SNPs 
(rs1041981, rs2239704, rs2229094 and rs746868) are present 
in high LD with rs909253 [27,28], inconsistent results of the 
four SNPs and cancer risk are observed for different cancers in 
Asian, North American and European populations 
[18-20,22-24,31-51]. In the current study, we perform a 
comprehensive meta-analysis to evaluate the effects of the four 
functional SNPs (rs1041981, rs2239704, rs2229094 and 
rs746868) on cancer risk. 

Materials and Methods 

Publication search 

The literatures included in the analysis were obtained from 
the databases of PubMed, Web of Science, and China National 
Knowledge Infrastructure (CNKI). The keywords applied the 
Medical Subject Headings (MeSH) in the US National Library of 
Medicine terms including "LTA", "Lymphotoxin alpha", "TNF- 
beta", "polymorphism", and "rs1041981", "Thr26Asn", 
"rs746868", "rs2229094", "Cys13Arg" or "rs2239704" paired 
with "cancer" or "tumor". Meanwhile, the manuscripts should be 
published in Chinese or English up to July 2013. Full texts were 
read to select the relevant information. The related articles in 
the MEDLINE as well as the reference lists of all retrieved 
studies were also checked for citations of other relevant 
publications that were not identified initially. 



Inclusion criteria 

Articles enrolled in our meta-analyses met the following 
inclusion criteria: (1) evaluating the association between LTA 
polymorphism rs1041981, rs2239704, rs2229094 or rs746868 
and cancer risk; (2) case-control study; (3) results with 
sufficient published data to estimate an OR with a 95% CI; (4) 
the genotype distribution in controls met Hardy-Weinberg 
equilibrium (HWE). 

Data extraction 

Two investigators (YH and XY) independently extracted the 
data from all eligible publications basis on the selection criteria 
listed above. Any disagreements were established by 
discussion until a consensus was reached. If there were 
numerous publications from the same study group, the most 
complete and recent results were extracted. Useful information 
collected from each study comprise the first author's surname, 
the year of publication, country, ethnic population, source of 
control groups (population- or hospital-based), cancer type, 
total sample size, genotyping method and the number of 
genotype distribution in cases and controls. 

Statistical analyses 

ORs and 99% CIs in the case-control studies were employed 
to assess the association between the LTA polymorphisms and 
cancer risk. The pooled ORs were performed under the 
additive, dominant, and recessive models. Z test was employed 
to estimate the significance of pooled OR. The departure of 
HWE for the genotype distribution in controls were analyzed by 
the Arlequin program (version 3.5) [53]. The meta-analyses 
were performed using the Review Manger and Stata software 
(version 11.0, Stata Corporation, College Station, TX) [54]. A 
chi-square based Q-statistic test was calculated for the 
heterogeneity of studies in the meta-analysis [55,56]. The 
inconsistency index (I 2 statistic) was examined to judge the 
heterogeneity between studies [56]. An I 2 > 50% suggested a 
significant heterogeneity in the meta-analysis. Random-effect 
[57] or fixed-effect [55,58] models were used for, the meta- 
analysis with significant heterogeneity (I 2 > 50%) or the one 
with minimal to moderate heterogeneity (I 2 < 50%), 
respectively. Subgroup meta-analyses were performed by 
cancer type, population, and source of control. A sensitivity 
analysis was also performed by excluding each study. The 
Begg's funnel plots and Egger regression test were used to 
evaluate the Publication bias [59]. A two-tailed P < 0.05 was 
considered statistically significant. 

Results 

Literature flow 

As shown in Figure 1 , the literature search identified a total 
of 248 potentially relevant records. After reading the title or 
abstract, 50 studies concerning the association for the four LTA 
polymorphisms and cancer risk were considered for the 
following step. A total of 16 articles were excluded for 
insufficient genotyping information, or ineligible samples, or 
other LTA polymorphism, or duplicated studies. Finally, we 
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248 potentially relevant records identified through database 
searching of PubMed. Web of Science and CNKI 



5 

202 publications assessed for eligibility 



-2 



151 studies unqualified: 
•Not polymorphism studies 
^ 'Not case-control studies 
•No! cancer studies 
•Reviews or letters 



5 1 studies concerning the association for 
LTA polymorphisms and cancer risk 
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27 studies excluded for: 
•Other polymorphism in L T.4 
•Removal of duplicates 
•Ineligible samples 


> 



24 articles (including 30 study stages) were finally included 
in meta-analysis: 

• 8 articles (including 1 1 stages) tot LTA rsI041981: 

• 10 articles (including 13 stages) for LTA rs2239704; 
5 articles (including 9 stages) for LTA rs222909-l: 

7 articles (including 7 stages) for LTA rs746868 



Figure 1. Flow diagram depicts literature search and 
study selection. 

doi: 10.1371/journal.pone.0082519.g001 

collected 24 articles [18-20,31-52] (including 30 study stages) 
focusing on the relationship of rs1041981, rs2239704, 
rs2229094 or rs746868 polymorphisms with cancer risk. 

Study characteristics 

Characteristics of enrolled studies were summarized in Table 
1. Among the 24 eligible studies, 8 articles 
[18,20,31,34,35,39,48,52] (including 11 stages with 7,483 
cases and 11,938 controls) focused on the association 
between rs1041981 and cancers. A total of 10 studies 
[19,33,36,38,40,41,43,46,47], reported the association of 
cancers with rs2239704 including 13 stages among 6,049 
cases and 7,621 controls. And 5 studies (including 9 stages) 
[19,31,46,51,52] among 7,133 cases and 10,305 controls 
evaluated the contribution of rs2229094 to the risk of cancers. 
There were 7 articles (including 3,487 cases and 4,633 
controls) [32,37,42,44-46,49] involved with the association 
between rs746868 and cancers. All the records were collected 
from English publications. There were 7 Asian studies, 9 
European studies, and 14 North American studies. Most of the 
cancer types were adenocarcinoma and hematological 
malignancy including gastric [31,34,42,45,46,49], breast 
[32,51,52], lung [31], prostate [38,44], and cervical cancers 
[35,39,48], colorectal adenoma [18,37], Non-Hodgkin 
Lymphoma (NHL) [19,20,33,36,40,47], Chronic Lymphocytic 
Leukaemia (CLL) [43] and Testicular Germ Cell Tumors 



(TGCT) [41]. The 30 selected studies in the meta-analysis 
included 13 studies with hospital-based controls, 17 studies 
with population-based controls. The controls were primarily 
population-based ones matched for ethnicity, age, gender or 
study region. 

As shown in Table 2, the genotype distribution of the four 
SNPs and the genotyping method of the collected studies were 
retrieved carefully. Genotype distributions of four 
polymorphisms in all the controls met HWE (P > 0.05). 

Meta-analysis of rs1041981 

The main results of the meta-analysis for rs1041981 
polymorphism were presented in Table 3. No heterogeneity 
existed in the involved studies (I 2 = 16.0%, P = 0.29). The 
summary effect OR was 1.10 for G allele (99% CI = 1.04-1.16, 
P = 0.001, Figure 2). Subgroup analysis suggested that 
rs1041981 increased the risk of several types of cancer, such 
as adenocarcinoma (OR = 1.07, 99% CI = 1.01-1.14, P = 
0.030, I 2 = 0.08%), squamous carcinoma (OR =1.19, 99% CI = 
1.03-1.37, P = 0.018, I 2 = 46.0%), hematological malignancy 
(OR = 1.73, 99% CI = 1.08-2.77, P = 0.023). The significant 
association were also found in multiple populations including 
Asians (OR = 1.18, 99% CI = 1.04-1.34, P = 0.010, I 2 = 0.0%) 
and Europeans (OR = 1.07, 99% CI = 1.00-1.14, P = 0.047, I 2 = 
0.0%). In addition, the source analysis showed positive 
association in both hospital-based group (OR = 1 .24, 99% CI = 
1.09-1.42, P = 0.001, I 2 = 0.0%) and population-based group 
(OR = 1.07, 99% CI = 1.00-1.14, P = 0.042, I 2 = 0.0%, Table 3). 

Meta-analysis result of rs2239704 

Figure 3 showed the forest plot for the rs2239704 stratified 
by cancer type. Significantly increased cancer risk was found in 
the G versus A model based on the studies (OR = 1.08, 99% 
CI = 1.01-1.16, P = 0.021, I 2 = 0.0%, Table 4). In the subgroup 
meta-analysis by cancer type, the rs2239704-G was observed 
with positive association with hematological malignancy (OR = 
1.10, 99% CI = 1.01-1.20, P = 0.023, I 2 = 0.0%, Table 4). 
Significant associations of rs2239704 with hematological 
malignancy cancers were also found in North American 
populations (OR = 1.09, 99% CI = 1 .01-1 .17, P = 0.025, I 2 = 
2.7%) and in the population-based group (OR = 1.08, 99% CI = 
1.00-1.17, P = 0.044, I 2 = 2.2%, Table 4). 

Meta-analysis result of rs2229094 

The main pooled data for rs2229094 polymorphism were 
listed in Table 5. For the overall data of the enrolled 9 study 
stages, significant association of rs2229094 polymorphism with 
cancer risk were shown on Figure 4 (OR = 1.28, 99% CI = 
1.09-1.50, P = 0.003, I 2 = 0.0%). In the additional analysis, 
significantly increased risks were observed in adenocarcinoma 
(OR = 1.33, 99% CI = 1.11-1.59, P = 0.002, I 2 = 0.0%). In the 
stratified analysis by population, strong association of 
rs2229094 with adenocarcinoma cancers was found in Asians 
(OR = 1.92, 99% CI = 1.04-3.57, P = 0.038, I 2 = 0.0%) and 
Europeans (OR = 1.26, 99%CI = 1.02-1.56, P = 0.029, I 2 = 
0.0%). Further subgroup analysis in recessive model, 
statistically associations were presented both hospital-based 
(OR = 1 .92, 99% CI = 1 .04-3.57, P = 0.038, I 2 = 0.0%) and 
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Table 1. Characteristics of studies included in the meta-analyses 3 . 



First author (year) 


Case/control Cancer type Populations 


SNP 


Source of control Matching 


Takei K (2008) 


570/581 


AC 


Asian 


rs2229094, rs1041981 


HB 


Ethnicity, study region 


Takei K (2008) 


159/581 


AC 


Asian 


rs2229094, rs1041981 


HB 


Ethnicity, study region 


Takei K (2008) 


177/581 


AC 


Asian 


rs2229094, rs1041981 


HB 


Ethnicity, study region 


Mahajan R (2008) 


301/415 


AC 


European 


rs2229094, rs2239704, rs746868 


PB 


Age, gender, study region, ethnicity 


Wang SS (2009) 


990/827 


HM 


North American 


rs2229094,rs2239704 


PB 


Ethnicity, frequency 


Wang SS (2009) 


434/515 


HM 


North American 


rs2229094,rs2239704 


PB 


Ethnicity, frequency 


Wang SS (2009) 


518/462 


HM 


North American 


rs2229094,rs2239704 


PB 


Ethnicity, frequency 


Abbas S (2010) 


3136/5477 


AC 


European 


rs2229094,rs1041981 


PB 


Age, study region, ethnicity 


Madeleine MM (2011) 


848/866 


AC 


North American 


rs2229094 


PB 


Age, ethnicity 


Wang SS (2006) 


1029/842 


HM 


North American 


rs2239704 


PB 


Age, gender 


Lan Q (2006) 


417/499 


HM 


North American 


rs2239704 


PB 


Age, ethnicity, residence 


Liu X (2006) 


822/861 


AC 


North American 


rs2239704 


PB 


Age, ethnicity, date of blood collection 


Purdue MP (2007) 


210/602 


OC b 


North American 


rs2239704 


PB 


Age, gender, ethnicity, date of collection 


Purdue MP (2007) 


296/602 


OC b 


North American 


rs2239704 


PB 


Age, gender, ethnicity, date of collection 


Purdue MP (2007) 


352/463 


HM 


European 


rs2239704 


HB 


Age, gender, residence 


Cerhan JR (2008) 


441/474 


HM 


North American 


rs2239704 


HB 


Age, gender, residence 


Jacobs EJ (2008) 


453/1184 


AC 


North American 


rs2239704, rs746868 


HB 


Age, ethnicity, date of blood collection 


Ennas MG (2008) 


38/112 


HM 


European 


rs2239704 


PB 


Age, gender 


Lee SG (2004) 


328/253 


AC 


Asian 


rs1041981 


PB 


Gender, study region 


Niwa Y (2005) 


44/320 


AC 


Asian 


rs1041981 


HB 


Gender, study region 


Niwa Y (2005) 


87/320 


SC b 


Asian 


rs1041981 


HB 


Gender, study region 


Niwa Y (2007) 


110/220 


AC 


Asian 


rs1041981 


HB 


Age, gender 


Aissani B (2009) 


135/140 


HM 


North American 


rs1041981 


HB 


CD4+ counts, the duration of HIV infection 


Castro FA (2009) 


951/1707 


SC b 


European 


rs1041981 


PB 


Age, ethnicity, study region 


Sainz J (2012) 


1760/1727 


AC 


European 


rs1041981 


PB 


Age, gender, study region 


Gunter MJ (2006) 


219/205 


AC 


North American 


rs746868 


HB 


Age, gender, ethnicity, study region 


Gaudet MM (2007) 


1622/1288 


AC 


North American 


rs746868 


PB 


Age, residence 


Garcia-Gonzalez (2007) 


404/404 


AC 


European 


rs746868 


HB 


Age, gender, ethnicity, residence 


Crusius JB (2008) 


428/1124 


AC 


European 


rs746868 


PB 


Age, gender, residence, date of collection 


Garcia-Gonzalez (2009) 


57/24 


AC 


European 


rs746868 


HB 


Ethnicity, study region 



a: PB, Population based; HB, Hospital based; AC: Adenocarcinoma; HM: Hematological malignancy; SC: Squamous carcinoma; OC: Other cancer 

b: The two studies of squamous carcinoma are all cervical cancer; the two studies of other cancer are testicular germ cell tumors (TGCT) including seminoma and no 
seminoma. 

doi: 10.1371/journal.pone.0082519.t001 



population based groups (OR = 1.24, 99% CI = 1.05-1.46, P = 5, P = 0.084 for rs1 041 981 ; P = 0.602 for rs2239704; P = 0.433 
0.01 0, I 2 = 0.0%, Table 5). for rs2229094; P = 0.343 for rs746868). 



Meta-analysis result of rs746868 

There was no evidence of association between rs746868 
and cancer risk (OR= 1.01, 99% CI = 0.93-1.10, P = 0.771, I 2 = 
0.0%, Table S1). Alternative genetic models and subgroup 
meta-analyses did not reveal any significant results (Table S1). 

Sensitivity analyses and publication bias 

The sensitivity analyses were performed by excluding each 
study. The statistical significance of the results was not 
changed (data not shown). The Egger's population bias plot 
indicated no visual publication bias in the meta-analysis (Figure 



Discussion 

In the systematic review and meta-analysis, we analyze 30 
case-control studies among 58,649 participants to verify the 
association between the four SNPs and cancer risk. Our 
findings suggest that rs1041981, rs2239704 and rs2229094 
are able to increase cancer risk by the significant association 
results under different genetic models. No evidence is found for 
the association between rs746868 and cancer risk. 

In subgroup meta-analyses stratified by cancer type, our 
results show that rs2239704 and rs2229094 are likely to 
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Table 2. Distribution of the four SNPs genotypes between cancer and control group included in the meta-analyses 0 . 



SNP 


First author (Year) 


Case (AA/Aa/aa) 


Control (AA/Aa/aa) 


HWE 


Genotyping 


rs2229094 


Takei K (2008) 


374/179/17 


396/172/13 


0.30 


Melting curve 


(TT/TC/CC) 


Takei K (2008) 


98/52/9 


396/172/13 


0.30 


Melting curve 




Takei K (2008) 


110/58/9 


396/172/13 


0.30 


Melting curve 




Mahajan R (2008) 


206/74/21 


247/150/18 


0.48 


TaqMan 




Wang SS (2009) 


530/387/73 


433/333/61 


0.86 


TaqMan 




Wang SS (2009) 


231/167/36 


269/209/37 


0.74 


TaqMan 




Wang SS (2009) 


282/179/39 


276/160/26 


0.69 


TaqMan 




Abbas S (2010) 


1686/1199/251 


2965/2153/359 


0.23 


MassArray 




Madeleine MM (2011) 


444/329/75 


475/334/57 


0.93 


SNPlex 


rs2239704 


Wang SS (2006) 


421/446/162 


332/370/140 


0.05 


TaqMan 


(CC/AC/AA) 


Lan Q (2006) 


165/189/63 


186/226/87 


0.23 


TaqMan 




Liu X (2006) 


285/537 


351/510 


> 0.05 


TaqMan 




Purdue MP (2007) 


82/99/29 


221/290/91 


0.86 


TaqMan 




Purdue MP (2007) 


89/154/53 


221/290/91 


0.86 


TaqMan 




Purdue MP (2007) 


138/168/46 


162/229/72 


0.63 


TaqMan 




Cerhan JR (2008) 


169/217/55 


170/225/79 


0.78 


TaqMan 




Ennas MG (2008) 


14/17/7 


36/53/23 


0.70 


TaqMan 




Mahajan R (2008) 


85/1 38/76 


105/223/85 


0.12 


TaqMan 




Wang SS (2009) 


357/373/132 


284/311/116 


0.06 


TaqMan 




Wang SS (2009) 


143/152/51 


162/189/74 


0.15 


TaqMan 




Wang SS (2009) 


197/229/58 


153/214/66 


0.62 


TaqMan 


rs1041981 


Lee SG (2004) 


109/156/63 


74/1 32/47 


0.44 


PCR-sequencing 


(CC/AC/AA) 


Niwa Y (2005) 


17/23/4 


107/165/48 


0.25 


AmpliTaq Gold 




Niwa Y (2005) 


43/36/8 


107/165/48 


0.25 


AmpliTaq Gold 




Niwa Y (2007) 


51/43/16 


71/114/35 


0.40 


AmpliTaq Gold 




Takei K (2008) 


211/280/87 


195/304/89 


0.11 


Melting curve 




Takei K (2008) 


73/71/20 


195/304/89 


0.11 


Melting curve 




Takei K (2008) 


74/73/32 


195/304/89 


0.11 


Melting curve 




Aissani B (2009) 


50/63/22 


71/61/8 


0.39 


BeadArray 




Castro FA (2009) 


341/456/154 


557/813/337 


0.20 


MassArray 




Abbas S (2010) 


1498/1317/332 


2481/2399/607 


0.47 


MassArray 




Sainz J (2012) 


833/729/198 


794/760/173 


0.69 


KASPar Assay 


rs746868 


Gunter MJ (2006) 


85/107/27 


76/102/27 


0.46 


TaqMan 


(GG/GC/CC) 


Gaudet MM (2007) 


576/769/277 


446/621/221 


0.86 


AmpFISTR 




Garcia-Gonzalez (2007) 


135/194/75 


142/191/71 


0.48 


TaqMan 




Mahajan R (2008) 


83/143/74 


108/220/84 


0.16 


TaqMan 




Crusius JB (2008) 


151/205/72 


398/545/181 


0.85 


TaqMan 




Garcia-Gonzalez (2009) 


19/26/12 


9/10/5 


0.67 


RFLP 



c: The results of Jacobs EJ et al (2008) are not listed in the table because they only proved the A allele frequency. The study of Liu X et al (2006) only proves the number of 
CC/CA+AA. 

doi: 10.1371/joumal.pone.0082519.t002 
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Table 3. Meta-analysis of rs1 041 981 . 





Variables 




Additive model 




Dominant model 




Recessive model 






N 


P, OR(99% CI) 


P(Q-test). I 2 


P, OR(99% CI) 


P(Q-test), I 2 


P, OR(99% CI) 


P(Q-test), I 2 


Total 


11 


0.001, 1.10(1.04-1.16) 


0.29, 16.0% 


<0.0001, 1.15(1.07-1.25) 


0.33, 12.2% 


0.035, 1.13(1.01-1.27) 


0.76, 0.0% 


Cancer type 


Adenocarcinoma 


8 


0.030, 1.07(1.01-1.14) 


0.69, 0.08% 


0.003, 1.14(1.04-1.24) 


0.43, 0.0% 


0.246, 1.08(0.95-1.23) 


0.99, 0.0% 


Squamous carcinoma 


2 


0.018, 1.19(1.03-1.37) 


0.17, 46.0% 


0.061, 1.21(0.99-1.49) 


0.12, 58.3% 


0.053, 1.30(1.00-1.69) 


0.56, 0.0% 


Hematological malignancy 


1 


0.023, 1.73(1.08-2.77) 


NA 


0.083, 1.75(0.93-3.29) 


NA 


0.040, 3.21(1.05-9.80) 


NA 


Population 


Asian 


7 


0.010, 1.18(1.04-1.34) 


0.64, 0.0% 


<0.0001, 1.39(1.16-1.66) 


0.35, 11.4% 


0.357, 1.12(0.88-1.43) 


0.93, 0.0% 


North American 


1 


0.023, 1.73(1.08-2.77) 


NA 


0.083, 1.75(0.93-3.29) 


NA 


0.040, 3.21(1.05-9.80) 


NA 


European 


3 


0.047, 1.07(1.00-1.14) 


0.38, 0.0% 


0.035, 1.10(1.01-1.20) 


0.84, 0.0% 


0.096, 1.12(0.98-1.28) 


0.52, 0.0% 


Source of control 


Hospital based 


7 


0.001, 1.24(1.09-1.42) 


0.44, 0.0% 


<0.0001, 1.45(1.20-1.74) 


0.68, 0.0% 


0.121, 1.23(0.95-1.60) 


0.61, 0.0% 


Population based 


4 


0.042, 1.07(1.00-1.14) 


0.58, 0.0% 


0.027, 1.10(1.01-1.20) 


0.92, 0.0% 


0.114, 1.11(0.98-1.26) 


0.66, 0.0% 



NA: not applicable for the heterogeneity test, 
doi: 10.1371/joumal.pone.0082519.t003 



Study 
ID 



rs1041981 Additive model - Cancer type 

OR (99%CI) %Weight 



Adenocarcinoma 

Lee SG et al (Korean-2004) 

Niwa Y et al (Japan-2005) 

Niwa Y et al (Japan-2007) 

Takei K et al (Japan-2008) 

Takei K et al (Japan-2008) 

Takei K et al (Japan-2008) 

Abbas S et al (Germany-2010) 

Sainz J et al (Germany-2012) 

Subtotal (l-squared = 0 0%, p = 0 694) 

Squamous carcinoma 
Niwa Y et al (Japan-2005) 
Castro FA et al (Sweden-2009) 
Subtotal (l-squared = 46.0%. p = 0.173) 

Hematological malignancy 
Aissani B et al (USA-2009) 
Subtotal (l-squared = .%. p = .) 

Heterogeneity between groups: p = 0.070 
Overall (l-squared = 16.0%. p = 0.292) 



o 



o 



1.07 (0 79,1 45) 
1.27 (0.69,2.33) 
1 39 (0.89.2 16) 
1.07 (0 86, 1.33) 
1.36 (0.97,1,90) 
1.12(0.81. 1.54) 
1.07 (0.98.1.16) 
1.00 (0.88,1.14) 
1.07(1.01,1.14) 



1.62 (1.01,2 59) 
1.15 (0.99, 1.33) 
1.19 (1.03.1.37) 



3 -46 

0.86 

1.62 

6.71 

2.82 

3.09 

44.87 

19.04 

82 48 



1.44 
14.64 
16 08 



1.73 (1.08,2.77) 1 44 
1.73(1.08,2 77) 1 44 



1.10 (1.04,1 16) 100.00 



Figure 2. Meta-analysis of LTA rs1041981 polymorphism and cancer risk in the additive model stratified by cancer 
type. 

doi: 10.1371/journal.pone.0082519.g002 



increase the risk of hematological malignancy and 
adenocarcinoma, respectively. Hematological malignancy is 
one cancer of uncontrolled proliferation of clonal B-cells, T-cell 
or NK-cell at different stages of maturations [60]. 
Adenocarcinoma can arise in many tissues of the body 
including stomach, breast, pancreas, colon and so on [61]. 
Both adenocarcinoma and hematological malignancy are multi- 
factorial diseases with complex interactions of genetic and 
environmental factors [62]. Accumulating evidences have 
reported that genetic variations in immune genes are 



susceptibility loci for hematological malignancy [47,63]. LTA 
rs2239704 polymorphism is the main predictor of LTA 
production in human B cells [28]. Adenocarcinoma is the most 
common histologic type of cancer. The LTA variants are shown 
to be associated with cancer risks under different genetic 
models. These genetic models consist of additive (per variant 
allele vs. common allele), dominant (variant allele carriers vs. 
homozygotes for the common allele) and recessive 
(homozygotes for the variant allele vs. all others) genetic 
models [64,65]. LTA rs1041981 is associated with a 
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Study 
ID 



rs2239704 Additive model - Cancer type 

OR (99%CI) %Weight 



Hematological malignancy 

Wang SS et al (USA-2006) 

Lan Q et al (USA-2006) 

Purdue MP et al (Australia-2007) 

Cerhan JR et al (USA-2008) 

Ennas MG et al (ltaly-2008) 

Wang SS et al (USA-2009) 

Wang SS et al (USA-2009) 

Wang SS et al (USA-2009) 

Subtotal (l-squared = 0.0%, p = 0 990) 

Adenocarcinoma 
Uu X et al (USA-2006) 
Jacobs EJ et al (USA-2008) 
Mahajan R et al (Poland-2008) 
Subtotal (l-squared = 55 .1%, p = 0 108) 

Other Cancer 

Purdue MP et al (USA-2007) 
Purdue MP et al (USA-2007) 
Subtotal (l-squared = 44.7%, p = 0.179) 

Heterogeneity between groups: p = 0 260 
Overall (l-squared = 0 0%, p = 0.598) 



o 



o 



1 05(0 88, 1 25) 
1.10(0 86, 1.41) 
1.15(0.88. 1.50) 
1.15(0.90. 1.47) 
1.15(0.57.2.30) 
1 05(0.87. 1.28) 
1.13(0 86. 1.49) 
1.20(0.94. 1 54) 
1.10(1.01. 1.20) 



14.42 

727 

6.25 

7.31 

091 

11.92 

5.88 

729 

61 25 



1.30(1 06.1.58) 11 15 

1 00(0 81.1.22) 1059 

0 96(0 73.1.27) 5 79 

1 10(0 97.1.25) 27.53 



1.08(0 80.1.46) 4 89 
0 82(0.63.1.07) 6 33 
0.92(0.76.1.13) 11 22 



1 08(1.01.1.16) 100.00 



Figure 3. Meta-analysis of LTA rs2239704 polymorphism and cancer risk stratified by cancer type. 

doi: 10.1 371 /joumal.pone.0082519.g003 



Table 4. Meta-analysis of rs2239704. 





Variables 




Additive model 




Dominant model 




Recessive model 






N 


P, OR(99% CI) 


P(Q-test). I 2 


P, OR(99% CI) 


P(Q-test), I 2 


P, OR(99% CI) 


P(Q-test), I 2 


Total 


13 


0.021, 1.08(1.01-1.16) 


0.60, 0.0% 


0.019, 1.11(1.02-1.21) 


0.76, 0.0% 


0.027, 1.12(1.01-1.24) 


0.70, 0.0% 


Cancer type 


Hematological malignancy 


8 


0.023, 1.10(1.01-1.20) 


0.99, 0.0% 


0.072, 1.12(0.99-1.26) 


0.99, 0.0% 


0.051, 1.17(1.00-1.38) 


0.99, 0.0% 


Adenocarcinoma 


3 


0.132 1.10(0.97-1.25) 


0.11, 55.1% 


0.036, 1.16(1.01-1.34) 


0.25, 27.6% 


0.135, 1.11(0.97-1.28) 


0.05, 65.8% 


Other cancer 


2 


0.440, 0.93(0.76-1.13) 


0.18, 44.7% 


0.414, 0.89(0.67-1.18) 


0.18, 45.1% 


0.691, 0.93(0.63-1.35) 


0.44, 0.0% 


Population 


North American 


10 


0.025, 1.09(1.01-1.17) 


0.41, 2.7% 


0.039, 1.10(1.00-1.21) 


0.53, 0.0% 


0.015, 1.14(1.03-1.27) 


0.75, 0.0% 


European 


3 


0.532, 1.06(0.88-1.28) 


0.64, 0.0% 


0.214, 1.19(0.90-1.57) 


0.99, 0.0% 


0.727, 0.94(0.68-1.31) 


0.38, 0.0% 


Source of control 


Hospital based 


3 


0.257, 1.08(0.94-1.24) 


0.61, 0.0% 


0.444, 1.07(0.90-1.26) 


0.72, 0.0% 


0.386, 1.09(0.90-1.31) 


0.44, 0.0% 


Population based 


10 


0.044, 1.08(1.00-1.17) 


0.42, 2.2% 


0.022, 1.13(1.02-1.24) 


0.60, 0.0% 


0.038, 1.14(1.01-1.28) 


0.61, 0.0% 



doi: 10.1371/journal.pone.0082519.t004 

significantly lower presence of Japanese male lung cancer 
under the dominant model (CA+AA versus CC) [28]; LTA 
rs909253 is associated with high risk of Asian gastric cancer in 
the heterozygote comparison (GA versus AA) [15]; LTA 
rs2009658 is shown to be associated with significantly elevated 
risk of breast cancer among Caucasian women aged 45-64 
years under the additive model (C versus G) [51]. Our results 
showed that LTA rs1041981 and rs2239704 polymorphisms 
were correlated with cancer in three genetic models, 
suggesting an additive effect for these two polymorphisms on 
the risk of cancer. In contrast, rs2229094 was associated with 
the risk of adenocarcinoma in recessive model, implying a lack 



of contribution for the heterozygote to the risk of 
adenocarcinoma. The results could partly be attributed to 
different LTA polymorphisms that play different roles in different 
cancers and different populations. Additionally, the different 
tumors in humans may be generated by special carcinogenic 
mechanisms which would lead to multiple connections with one 
same genetic locus [66]. 

Our subgroup meta-analyses by ethnicity find that rs1041981 
shows a negative relationship with cancers in Asians under the 
recessive model. SNP rs2239704 shows positive association in 
North Americans but not Europeans. Additionally, rs2229094 is 
associated with cancer in Asians but not in North Americans or 
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Table 5. Meta-analysis of rs2229094. 





Variables 




Additive model 




Dominant model 




Recessive model 






N 


P, OR(99% CI) 


P(Q-test). I 2 


P, OR(99% CI) 


P(Q-test). I 2 


P, OR(99% CI) 


P(Q-test). I 2 


Total 


9 


0.076, 1.06(0.99-1.13) 


0.55, 0.0% 


0.476, 1.03(0.95-1.12) 


0.46, 0.0% 


0.003, 1.28(1.09-1.50) 


0.82, 0.0% 


Cancer type 


Adenocarcinoma 


6 


0.064, 1.07(1.00-1.15) 


0.37, 7.0% 


0.415, 1.04(0.95-1.14) 


0.25, 24.8% 


0.002, 1.33(1.11-1.59) 


0.70, 0.0% 


Hematological malignancy 


3 


0.762, 1.06(0.99-1.13) 


0.57, 0.0% 


0.996, 1.00(0.85-1.18) 


0.64, 0.0% 


0.456, 1.13(0.82-1.57) 


0.70, 0.0% 


Population 


Asian 


3 


0.041, 1.23(1.01-1.50) 


0.65, 0.0% 


0.104, 1.21(0.96-1.52) 


0.79 0.0% 


0.038, 1.92(1.04-3.57) 


0.82, 0.0% 


European 


2 


0.457, 1.03(0.95-1.13) 


0.19, 42.4% 


0.868, 0.99(0.89-1.11) 


0.06, 71.4% 


0.029, 1.26(1.02-1.56) 


0.52, 0.0% 


North American 


4 


0.347, 1.06(0.94-1.18) 


0.62, 0.0% 


0.645, 1.03(0.90-1.19) 


0.71 , 0.0% 


0.168, 1.21(0.92-1.58) 


0.76, 0.0% 


Source of control 


Hospital based 


3 


0.041, 1.23(1.01-1.50) 


0.65, 0.0% 


0.104, 1.21(0.96-1.52) 


0.79 0.0% 


0.038, 1.92(1.04-3.57) 


0.63, 0.0% 


Population based 


6 


0.244, 1.04(0.97-1.12) 


0.61, 0.0% 


0.877, 1.01(0.92-1.10) 


0.40, 1 .9% 


0.010, 1.24(1.05-1.46) 


0.90, 0.0% 



doi: 10.1371/joumal.pone.0082519.t005 



Study 
ID 



rs2229094 Recessive model - Cancer type 



Adenocarcinoma 
Takei K et al (Japan-2008) 
Takei K et al (Japan-2008) 
Takei K et al (Japan-2008) 
Mahajan R et al (Poland-2008) 
Abbas S et al (Germany-201 0) 
Madeleine MM et al (USA-201 1) 
Subtotal (l-squared = 0.0%. p = 0.703) 

Hematological malignancy 
Wang SS et al (USA-2009) 
Wang SS et al (USA-2009) 
Wang SS et al (USA-2009) 
Subtotal (l-squared = 0 0%. p = 0.698) 

Heterogeneity between groups: p = 0 406 

Overall (l-squared = 0.0%. p = 0.820) 



1 22 



o 



I 



OR (99%CI) %Weight 



1 3-4 (0.51 .3.51) 
2.62 (0.84,8.21) 
2.34(0.75,7.31) 
1.65 (0.71,3.87) 
1.24(1.00, 1.54) 
1 38(0.86.2.20) 
1.33(1.11, 1.59) 



1 00(0 63.1.59) 
1.17(0.62.2.18) 
1.42 (0.72,2.78) 
1.13(0.82. 1.57) 



8 21 



2 73 
1.95 

I. 98 
3.53 
54.47 

II. 51 
76.16 



11.85 
6.42 
5.56 
23.84 



128(1.09,1.50) 100.00 



Figure 4. Meta-analysis of LTA rs2229094 polymorphism and cancer risk stratified by cancer type. 

doi: 10.1371/journal.pone.0082519.g004 



Europeans. This might be explained by the different genotype 
and allele frequencies of these SNPs in subjects with different 
clinical characteristics, geographic distributions, and ethnic 
descent. Therefore, we could not exclude that the negative 
association in Europeans or North Americans may be due to a 
lack of power. Thus, matching criteria and selection bias and 
the stages of the cancers should be considered in the future 
case-control studies. Large sample size and various 
populations study would get more believable result in the 
future. 

In the present study, we have collected a series of 
parameters (including cancer type, population, source of 
control, genotyping method, and matching condition) to yield 



reliable result in the meta-analyses. We require that the 
genotype distribution in the controls met HWE (P > 0.05). We 
also perform subgroup meta-analyses by the collected 
parameters to reduce the potential stratification among the 
involved case-control studies. However, there are some 
limitations in our meta-analysis should be mentioned. Firstly, 
our meta-analyses combine the genetic studies from various 
cancers that may introduce dramatic stratifications in the meta- 
analyses, although we have controlled several parameters. 
Secondly, Our results can't tell the true causal variant of cancer 
risk from the three significant SNPs (rs1041981, rs2239704, 
rs2229094) since they are in high LD. Moreover, it is possible 
that the true causal variant may be other variant in high LD with 
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Figure 5. Funnel plots of LTA polymorphisms and cancer risk data for publication bias. 

doi: 10.1371/joumal.pone.0082519.g005 



them. Future investigation for the mechanistic roles of them is 
needed. Thirdly, there are only a few studies in the African 
populations. Future study in African populations needed to be 
performed in a large size cohort to investigate whether the 
negative results of subgroup meta-analysis in the African 
populations is due to a lack of power or genetic heterogeneity. 

In summary, the overall data of the present analyses suggest 
that three LTA variants (rs1041981, rs2239704 and rs2229094) 
can significantly increase the risk of cancers. Further well- 
designed studies in view of these variants are needed to 
explore their mechanistic roles in the pathogenesis of cancers. 
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